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CN ■ We investigate the condensation of charged K meson and neutral K° meson 

in dense neutron star matter. Calculations are performed in relativistic mean 
field models in which both the baryon-baryon and (anti)kaon-baryon interac- 
tions are mediated by meson exchange. It is found that K° condensation is 
quite sensitive to the antikaon optical potential and depends more strongly 
C\) ■ on the nucleonic equation of state. For moderate values of antikaon potential 

, and a rather stiff equation of state, a significant region of maximum mass 

ff~j \ star will contain K° meson. The critical density of K° condensation is always 

higher than that of K~ condensation. With the appearance of K~ and K° 
q ■ condensates, pairs of p — K~ and n — K° are produced with equal proportion 

■ leading to a perfectly symmetric matter of nucleons and antikaons in neutron 

stars. Along with K~ condensate, K° condensate makes the equation of state 
much softer resulting in smaller maximum mass stars compared to the case 
without any condensate. 
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I. INTRODUCTION 

The study of dense matter in laboratories reveals many new and interesting results. There 
is a growing interplay between the physics of dense matter and the physics of compact objects 
The composition and structure of a neutron star primarily depends on the nature of 
strong interaction. The equation of state of neutron star matter encompasses a wide range 
of densities, from the density of iron nucleus at the star's surface to several times the normal 
nuclear matter density encountered in the core. Since the chemical potentials of nucleons 
and leptons increase rapidly with density in the stars interior, several novel phases with 
large strangeness fraction may appear. Among these possibilities, hyperon formation is 
quite a robust mechanism. This is possible when the neutron chemical potential becomes 
sufficiently large so that the neutrons at their Fermi surfaces can decay into hyperons via 
weak strangeness non- conserving interactions. The threshold densities for hyperon formation 
strongly depend on the nuclear equation of state 

In recent years considerable interest has grown in the study of the properties of kaons, 
K, and antikaons, K, in dense nuclear matter and in neutron star matter. Within a chi- 
ral SU(3)l x SU(3)r Lagrangian, it was demonstrated by Kaplan and Nelson || that a 
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negatively charged antikaon K~ may undergo Bose-Einstein condensation in dense baryonic 
matter formed in heavy ion collisions. It was further predicted in the chiral perturbation 
theory that a K~ condensed phase may form in the core of neutron stars 0-§| in consonance 
with kaon-nucleon scattering data |§ and K~ atomic data ||10|| . In these chiral Lagrangians, 
the (anti)kaons are directly coupled to nucleons. The strongly attractive i^~-nucleon inter- 
action increases with density lowering the effective mass m* K of the antikaons. Consequently, 
the in-medium energy of K~ meson, u>k-, decreases and s-wave K~ condensation sets in 
when uj k - equals to the K~ chemical potential \ik- which, in turn, is equal to the electron 
chemical potential p e for a cold catalyzed (neutrino-free) neutron star matter. The typical 
critical density for K~ condensation in nucleons-only star matter is about (3— 4)n , where n 
is the normal nuclear matter density. The exact value, however, depends on the model used, 
i.e. on the nucleonic equation of state and on the parameters employed, especially on the 
depth of the attractive K~ optical potential. The net effect of K~ condensation in neutron 
star matter is that K~ replaces electrons in maintaining charge neutrality and the energy of 
the condensate is lowered because of strongly attractive interactions between nucleons and 
K~ mesons. Due to the softening of the equation of state the masses of the stars are reduced 
in presence of K~ condensate |[7|jrT|--[14| . It was also found that in presence of hyperons, K ~ 
condensation was delayed to higher density, and may not even exist in maximum mass stars. 
Protoneutron (newly born and neutrino rich) stars with K~ condensate was shown to have 
maximum masses larger than those of cold catalyzed (old and neutrino free) stars - a reversal 
from ordinary nucleons-only star. These stars may undergo supernovae explosion and then 
could collapse to small mass black holes during deleptonization |0J. 

The situation for kaons is however quite different. Theoretical investigations based on 
Nambu— Jona-Lasinio model | 16fl ) chiral perturbation theory ||, and one-boson-exchange 



model | |T2l , |T3"|| yield a repulsive optical potential for K + in nuclear medium. Therefore, K + 
condensation is expected to be forbidden in neutron stars. 

In the traditional meson-exchange picture |T7[, (anti)kaons and baryons interact via the 
exchange of a, lu, and p mesons. In addition to the decrease of m* K with density here, the 
energy ujk- is lowered since K~ meson experiences attractive vector u;-meson potential due 
to G-parity of antikaons. In contrast, the isovector p-meson field is repulsive for K~ and 
inhibits K~ condensation delaying it to higher densities. In neutron stars at densities above 
normal nuclear matter values, most of the protons convert into neutrons by inverse /3-decay 
since it is energetically favorable to posses a small number of electrons. The ratio of electrons 
to nucleons which is again equal to the proton fraction due to charge neutrality is typically 
Ye — Y p « 0.1 — 0.15. Therefore, neutron star matter is highly asymmetric and quite distinct 
from ordinary symmetric nuclear matter. This suggests that the p-meson field should have 
a strong repulsive effect on K~ meson. On the other hand, the p-meson field induces an 
attractive field for K° meson which is an isodoublet partner of K~ meson. This should lower 
the effective energy of K° meson, u^o, compared to that of K~ meson thereby making K° 
meson condensation more favorable in neutron star matter. The critical density for s-wave 
neutral K° meson condensation is governed by the condition u^o = 0, and should depend 
sensitively on the equation of state (EOS). Using variational chain summation method with 



two-nucleon and three-nucleon interactions, Akmal et al. |18] recently predicted a neutral 
pion condensation phase transition in neutron star matter at a density of ~ 0.2 fm -3 . 
So far no calculation of neutral K° condensation and its impact on the gross properties 
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of neutron stars has been performed. In this paper, we investigate the effect of antikaon 
condensation with more emphasize on K° condensate on the constitution and structure of 
neutron star matter in the standard meson exchange model. For this purpose, we employ 
the usual relativistic mean field Lagrangian JlT]] for baryons interacting via meson exchanges 
and includes also the self-interaction of the mesons. The antikaons are treated in the same 
footing as the baryons which interact by the exchange of the same mesons. For the K~ 
meson, the Lagrangian density of Ref. [|TJ]] is used, and it is extended to include the K° 
meson. We shall demonstrate within this model that, apart from the K~ meson (as already 
studied in Ref. fT^TJ), the K° condensate may also exist inside a neutron star and has a 
significant influence on the star properties. We shall also show that the threshold densities 
for the antikaon condensates are quite sensitive to the poorly known high density regime of 
the equation of state. 

The paper is organized as follows. In section II we describe the relativistic mean field 
(RMF) model of strong interactions. The relevant equations for neutron star matter with 
antikaon condensate are summarized in this model. In section III the parameters of the 
model are discussed and effects of antikaon condensates in neutron star matter are presented. 
Section IV is devoted to summary and conclusions. 



II. THE FORMALISM 

The starting point in the present approach is a relativistic field theoretical model of 
baryons and (anti)kaons interacting by the exchange of scalar a, isoscalar vector u, and 
vector isovector p mesons. The total hadronic Lagrangian can be written as the sum of the 
baryonic, kaonic and leptonic parts, i.e. £ = Cb + Ck + ^-i- Considering all the charge states 
of the baryon octet B = {n,p, A, S + , E~, S°, E~, 2°}, the baryonic Lagrangian is given by 

C B =Y.^B (il^ -m B + g aB Cr - gu>Bl^ ~ \flpBlpjB ■ p^j IpB 

+~ (d^a - mla 2 ) - U{a) 

+ \mlu^ - ip^ • pT + \m% Pil ■ p» . (1) 

Here ipB denotes the Dirac spinor for baryon B with vacuum mass m B and isospin operator 
Tb- The scalar self- interaction term 

U(a) = \ 9 2^ + , (2) 



is included 0] to achieve a realistic compression modulus at normal nuclear matter density. 



The Lagrangian density for (anti)kaons in the minimal coupling scheme is given by [T^ 

C K = D^KD^K - m%KK , (3) 

with a covariant derivative = <9 M + ig^K^p, + ig P KfK • P M - The isospin doublet for the 
kaons is denoted by K = (K + ,K°) and that for the antikaons is K = (K~,K°). The 
effective mass of (anti)kaons in this minimal coupling scheme is given by 
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m* K = m K - g aK a , (4) 

where mx is bare kaon mass. In the mean field approximation (MFA) adopted here, the 
meson fields are replaced by their mean values and the baryon currents by those generated 
in the presence of the mean meson fields. For a uniform and static matter within MFA, only 
the time-like components of the vector fields, and the isospin 3-component of p-meson field 
have non-vanishing values. The mean meson fields are denoted by a, ujq, and po3- 

For s-wave (k = 0) condensation of antikaons K, the dispersion relation representing the 
in-medium energies of K = (K~, K°) is given by 

w K -, K° = m K 9uK^0 =F ^9pKp03 , (5) 

where the isospin projection I 3R = =Fl/2 for the mesons K~ (— sign) and K° (+ sign) are 
explicitly written in the expression. Since the a and uj fields generally increase with density 
and both being attractive for antikaons, the effective energies of K are lowered in nuclear 
medium. Moreover, in nucleons-only matter p 3 = n p ~ n n { n p and n n are the proton and 
neutron densities) is negative, thus the p-meson field inhibits K~ condensation whereas it 
favors K° condensation. Employing G-parity which simply transforms the sign of vector 
potential, the energies of kaons K = (K + , K°) are expressed as 

ko = m* K + g^R^o ± ~g P KPo3 ■ (6) 

The tu-meson field being repulsive for the kaons and dominates over the attractive u-meson 
field at high densities is suggestive of the fact that kaon condensation should be highly 
restricted. 

The meson field equations in presence of baryons and antikaon condensates can be derived 
from Eqs. (l)-(3) as 

dU 

mla = -— + Y.9*Bn s B + g aK J2 n K i (7) 

B R 

m 5^0 = 9u)B n B ~ 9uK J2 H K ' ( 8 ) 
B K 



m 



P03 = Y 9 P BhBn B + g P K Y hit^R ■ (9) 



P' 

B 



Here the scalar and number density of baryon B are respectively 

Q 2J B + 1 f k F B m 



n 



2vr 2 Jo (P + mf)V2 



5^fc 2 dk , (10) 
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n B = (2J B + l)^ , (11) 

with effective baryonic mass m* B = m B — g a Bcr, Fermi momentum fcp s , spin J B , and isospin 
projection I% B . Note that for s-wave K condensation, the scalar and vector densities of 
antikaons are same, and in the MFA is given by f!4 



n K -,Ro = 2 ( uj K - tR o + gu K u ± ^g P xPo3 ) KK = 2m* K KK . { 12 ) 
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The total energy density e = Eb+sr nas contributions from baryons, leptons, and antikaons. 
The baryonic plus leptonic energy density is 

e B = \ m W + \d2(y Z + ^s^ 4 + \rnlul + ^m 2 p p 2 03 

+ E^^ l kFB (k 2 + mfflVdk + Y.- 2 [ kl (k* + m^kUk. (13) 

■q ATI JO ^ 7T «/ 

The last term corresponds to the energy of leptons as required in a neutron star matter. 
The energy density for antikaons is 

£R = m *K ( n K- + n R0 ) . (14) 

Since antikaons form s-wave Bose condensates, they do not directly contribute to the pressure 
so that the pressure is due to baryons and leptons only 



I99 S 4 ^22 2 2 

-2 m * a --^o + 2 m ^o + 2^03 

+ 3 \ 2tt 2 Jo (P + mf)V2 + 32^^y o ( fc 2 + m 2)i/2 



(15) 



The pressure due to antikaons is contained entirely in the meson fields via their field equa- 
tions (7)-(9). 

At the interior of neutron stars, nucleons and electrons undergo usual /3-decay processes 
n — > p + e~ + v e and p + e~ —> n + v e . When the electron chemical potential becomes equal 
to the muon mass, electrons are converted to muons by e~ — > /i _ + + z/ e . In the present 
study of cold neutron star we may assume that the neutrinos have left the system freely. 
Therefore the chemical potentials of nucleons and leptons are governed by the equilibrium 
conditions 

A*n - Hp = He = Hp , (16) 

where /i n , /i p , /i e , and /i^ are the chemical potentials of neutrons, protons, electrons, and 
muons with /i n p = [k 2 Fn + m^) 1 / 2 + g^N^o + hN9 P NPo3- With the onset of K condensation, 
the strangeness changing processes that may occur are N ^ N + K and e~ # K~ + z/ e , 
where N = (n,p) and K = (if - , -ft' ) denote the isospin doublets for nucleons and antikaons, 
respectively. The requirement of chemical equilibrium yields 

/in - /ip = UK- = /ie , (17) 
^0 = 0, (18) 

where hk- and /j, k o are respectively the chemical potentials of K~ and K°. The above 
conditions dictate the onset of antikaon condensations. When the effective energy of K~ 
meson, lvk-, equals to its chemical potential, Hk~, which in turn is equal to the electro- 
chemical potential /z e , a K~ condensate is formed. While K° condensation follows when 
its in-medium energy satisfies the condition lvro = fi K o = 0. When other baryons in the 
form of hyperons are present in the neutron star matter, the standard /5-decay processes for 
nucleons generalize to the form B x — > B 2 + I + vi and B 2 + I — > B\ + ui, where B\ and B 2 
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are baryons and I is a lepton. All the equilibrium conditions involving the baryon octet may 
then be summarized by a single generic equation 



= ~ QBf^e , (19) 

where \xb and qs are, respectively, the chemical potential and electric charge of baryon 
species B. The above relations indicate that, in chemical equilibrium, two independent 
chemical potentials, \i n and /i e , exist which correspond to baryon number and electric charge 
conservation. For neutron star matter we need to include also the charge neutrality condi- 
tion, which in presence of antikaon condensate is expressed as 

VBriB - n K - - n e - = . (20) 

B 

In the present calculation we consider antikaon condensation as a second order phase 
transition. In principle in neutron star matter with two conserved charges, baryon and 
electric charge, the condensation can also be treated as a first order phase transition |14| . 



In this situation, the star would have at low density a normal phase of baryons and leptons 
followed by a mixed phase of K~ condensate and baryons, and possibly a pure K~ condensed 
phase at a higher density. However in the present situation this would be more complicated 
as K~ and K° have to be treated as two separate phases, and especially when the threshold 
density of one condensate may be reached in the mixed phase of nucleons and the other 
condensate. Treatment of such a triple phase would be numerically very involved and is 
postponed to a future publication. 



III. RESULTS AND DISCUSSIONS 

In the effective field theoretic approach adopted here, three distinct sets of coupling 
constants for nucleons, kaons, and hyperons associated with the exchange of a, u, and p 
mesons are required. The nucleon-meson coupling constants generated by reproducing the 
nuclear matter saturation properties are taken from Glendenning and Moszkowski of Ref. 
pi| . This set is referred to as GM1 and listed in Table I. 

Let us now determine the kaon-meson coupling constants. According to the quark and 
isospin counting rule, the vector coupling constants are given by 

9u>k = ^9luN and g pK = g pN . (21) 

The scalar coupling constant is obtained from the real part of the K~ optical potential at 
normal nuclear matter density 

Ur (no) = -9aK<7 - 9uK^0 ■ (22) 

The negative sign in the vector meson potential is due to G-parity. The critical density of 
K condensation should therefore strongly depend on the K~ optical potential. Recent fits 
P| f2T|j2^ | to the K~ atomic data indicate a strong attractive potential in dense matter. More 
recently, in high energy heavy ion collisions enhanced subthreshold K~ production with 
a steep spectral slope was found ||23| . In a simplified model this may be interpreted as a 
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strong K~ attraction in nuclear matter. On the other hand, available K~N scattering length 
suggests a repulsive interaction. The reason for this apparent ambiguity is the presence of 
A(1405)-resonance which is considered to be an unstable KN bound state just below the K~p 
threshold that makes the interpretation of the data more complicated. Recent analysis |22| 



of K~ atomic data using a hybrid model which combines the relativistic mean field approach 
in the nuclear interior and a phenomenological density dependent potential at low density 
showed that the real part of K~ optical potential could be as large as U K = —180 ±20 MeV 
at normal nuclear matter density while being slightly repulsive at low density in accordance 
with the low density theorem. Coupled channel formalism which automatically generates the 
A(1405)-resonance and successfully describes the low energy K~p scattering data yields an 



attractive potential for the K~ meson of U K (no) ~ —100 MeV |24]. In a chirally motivated 
coupled channel approach, the K optical potential depth was found to be U K (no) = —120 
MeV [^5| . The reason for such a spread in the predicted values lies in the diversity of treating 
the A(1405)-resonance. We have therefore determined the K — o coupling constant g a K for 
a set of values of U K (no) starting from —100 MeV to —180 MeV. This is listed in Table II 
for the set GM1. Since the u;- meson potential for K in this model is Kf (no) = g^K^a ~ 72 
MeV, a rather large sigma-kaon coupling constant of g a K = 3.674 is required to reproduce 
a depth of —180 MeV. Note that only for this large depth the value of the scalar coupling 
is similar to the prediction in the simple quark model i.e., g a x = gaN/3. 

As an alternative approach the kaon-meson coupling constants may be determined from 
the s-wave kaon-nucleon (KN) scattering length. At the tree level, the isospin averaged 
KN scattering length is given by ^6 



1 7=0 , 3 7=1 m K ( gaI<gaN gWguV \ fo „s 

a-KN - -7% + t^o - A n i 7 V 2 1 2 — • ^ Z6 

4 4 47r(l + mK/TriMj \ mi rnt J 

From the low density theorem, the kaon optical potential depth is given by 

27r / nw\ , . 

U K = (l + — )a KN n B . 24 

m K V m N J 

Using the experimental values [27| of Oq =1 = —0.31 fm and a^ = —0.09 fm, the scattering 



length is cirn — —0.255 fm. The K optical potential depth is then Uk ~ +29 MeV at 
n = 0.153 fm -3 and is repulsive. In the GM1 parameter set, for K depths of U K (no) = 
— 100(— 180) MeV, we obtain KN average scattering lengths and K optical potential depths 
of a KN = -0.468(-0.319) fm and U K (n ) ~ +44(-36) MeV. With U K (n ) = -180 MeV, 
the attractive potential obtained for kaon optical potential depth (apart from rather small 
vector potential of (anti)kaon V^f ) is due to the neglect of A(1405)-resonance which is beyond 
the scope of this paper. 

In the present study we shall be mostly concerned with antikaon condensation in 
nucleons-only matter. In general, the presence of hyperons delays the onset of K con- 



densation to much higher densities P,p!T|-4T3[] . However, for orientation we shall only provide 
briefly the effects of hyperons on antikaon condensation in neutron star matter. The vector 
coupling constants for the hyperons are obtained from SU (6) symmetry 

1 1 1 

T^giuN — 2#">A — — g^S , 

9 p n = = gpE ; g P A = . (25) 
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The cr-meson couplings to the hyperons (Y) can be obtained from the well-depth of Y in 
saturated nuclear matter 



U Y (no) = -Qayo- + g U Y^o • (26) 



Analysis of energy levels in A-hypernuclei suggests |28[ a well-depth of A in nuclear matter 
of Uj^(no) ~ —30 MeV. For the £ potential the situation is unclear, since there is no 



evidence for bound S-hypernuclei. The prediction range from completely bound E's |29] with 
Ug (n ) ~ -30 MeV to unbound gTJ with Ug (n ) ~ +30 MeV. A few events in emulsion 
experiments with K~ beams have been attributed [^] to the formation of H~-hypernuclei 
from which the depth of H in nuclear matter has been extracted to be Ug(n ) w -28 MeV. 
Recently, a few S-hypernuclei events have been identified in (K~,K + ) reaction from which 
a tighter constraint on S well-depth of ~ — 18 MeV has been imposed |5D , 3"I ]. The coupling 
constants g a y for hyperons have been obtained using these depths. 

We now present results for neutron star matter containing nucleons, leptons, and K 
condensates for the parameter set GM1 of Tables I and II. In Fig. 1, the nucleon scalar 
and vector potentials, and the electron chemical potential are displayed as a function of 
baryon density normalized to the equilibrium value of no for K~ optical potential depth 
of Ux{no) = —120 MeV. The solid lines refer to results for nucleons-only matter and the 
dashed ones correspond to calculation with K condensate. With K condensation we find 
that the scalar potential is enhanced while vector u- and p-meson potentials are decreased. 
The kinks in the dashed lines correspond to the onset of K° condensation when the meson 
potentials deviate further from nucleons-only results. The variation in the meson fields may 
be attributed to the occurrence of the source terms due to K condensation in the field 
equations of motion (Eqs. (7)-(9)). 

The populations of neutron star matter with and without K condensation are shown in 
Fig. 2. In the top panel, the particle abundances of nucleons-only matter are shown. Here all 
the particle fractions increase with baryon density. The behavior of the proton (and lepton) 
fractions are determined by nuclear symmetry energy which in turn is controlled by the 
p-meson field in the RMF models. In the central panel we exhibit the particle abundances 
of the star matter with K~ condensation at Ugino) = —120 MeV, which in fact is the first 
particle to condense among the two antikaons K = (K~,K°). Once K~ condensate sets 
in at 3.05n , it rapidly increases with density replacing the leptons in maintaining charge 
neutrality. The proton density which then becomes equal to K~ density eventually turns out 
to be larger than the neutron density; the neutron density nearly freezes to a constant value. 
This can be understood from the threshold condition of K~ condensation uo^- = p« — p P 
(see Eq. (17)). Substituting explicitly these values, we obtain 



(k 2 Fp + m^y /2 = (k 2 Fn + m*^y /2 + (-m* K + -g uN uJo ~ ^QpNPos) ■ (27) 



,V2 , , 9 aV2 , ( * , 1 1 

Since at high densities, the cu-meson potential dominates the cx-meson potential and 
P03 = rip — n n is negative, the term within the bracket is positive. It is therefore evi- 



dent that k,F v > kp„- This has been also observed in Ref. |I4J. Hence the net effect of 
K~ condensation is to lower the electron chemical potential p e and the nuclear symmetry 
energy (see also Fig. 1). In the bottom panel of Fig. 2, we show the particle fractions with 
both K~ and K° condensates. Interestingly, with the onset of K° condensate at 4.83no, the 
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neutron and proton abundances become identical resulting in an isospin saturated symmet- 
ric nuclear matter. This can be easily understood by substituting the threshold condition 
for K° condensation, i.e. u^o = 0, from Eq. (5) into Eq. (27). From these relations it is 
clear that K° condensate formation enforces the condition k Fp = k Fn . Alternatively, one 
may interpret from the central panel of Fig. 2 that once the neutron and proton densities 
(Fermi momenta) become identical, the isodoublet partner of K~ meson, i.e. K° should be 
populated. The electrochemical potential \x e = p„ — /j, p then merges with the magnitude of 
the p-meson potential \V^\ = |<7 p jvPo3| with p 03 = n^o — nx-, as is evident from Fig. 1. At 
densities above 7.32n , the antikaon densities are equal which is a general feature of isospin 
driving force towards symmetry. The p-meson field and thereby p e vanishes. As a conse- 
quence we are left with a system composed of exactly equal amount of n — p — K~ — K°, i.e. 
a perfectly symmetric matter of nucleons and antikaons. For this matter, the strangeness 
fraction f s = \S\/B = {jir- + n R°)/ n B is as large as unity. 

In Fig. 3, the effective masses of the (anti)kaons, m* K , are shown as a function of 
normalized density for various values of K~ optical potential depths from —100 MeV to 
— 180 MeV in steps of 20 MeV. Note that K is only a test particle in the field of nucleons 
and appear physically only after condensation. The decrease of m* K is found to be quite 
sensitive to the K~ depth which determines the onset of condensation. At a given Uk(tio), 
the (anti)kaon mass decreases more strongly with successive condensation of K~ and K° 
mesons. This is a manifestation of increase in the cx-meson field with condensation as can 
be seen in Fig. 1. 

Figure 4 shows the s-wave condensation energies for K~ meson, uk- (solid lines), and for 
K° meson, u^o (dashed lines), as a function of density of various antikaon optical potential 
depths in the GM1 set. The effective energies follow a similar qualitative trend as m* K . 
However, for depths up to U^fao) = —140 MeV, the effective kaon masses do not vary 
significantly (see Fig. 3). Therefore, considering Eq. (5), the large drop in lj k - xo with 
density is dominated by cu-meson that influences K condensation. For higher values of 
Uk, the large decrease in m* K s due to large couplings g a x are primarily responsible for 
condensation. For any value of U^, the isovector potential of nucleons shifts the energy 
of K° meson below that of the K~ meson by ~ 120 MeV. At a given U^, condensation 
of K~ meson occurs when u>k- intersects the electron chemical potential p e in absence of 
a condensate (shown by dotted line). For all values of K depth we find that K~ meson 
is formed before the threshold density for K° condensate, i.e. u^o = 0, is reached. The 
critical densities for K condensates are collected in Table III for different values of J7^(n ). 
A careful investigation of Fig. 4 reveals that beyond the densities for K~ condensates (and 
before K° threshold densities), the rate of decrease of the energies ou K - are attenuated. 
This can be traced back to the density-dependent behavior of the meson fields in Fig. 1. 
Above K~ formation density the decrease in tu-meson field dominates over the increase in 
the cr-meson field. This leads to smaller rate of decrease of ujk- and cu^o. Moreover, the 
decrease in the p-meson field due to K~ condensation causes even a smaller decrease of K° 
energy with density. Thus the K~ meson condensate with the help of the vector mesons 
conspires to delay the appearance of its isospin partner K° meson to a higher density. At 
densities above K° condensate, the K~ energy again starts to decrease due to dramatic fall 
of p-meson potential. 

The equation of state (EOS), pressure P versus the energy density e is displayed in Fig. 
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5 for the GM1 parameter set. The solid line represents calculation for nucleons-only star 
matter while the dashed lines refer to those with K formation for different values of U R (no). 
The strong attraction imparted by antikaon condensation makes the EOS softer. The kinks 
at higher densities in the dashed lines correspond to the appearance of K° condensation 
which further soften the EOS. Also the softness is quite sensitive to the choice of the K 
optical potential depth. A large U R (n ) corresponds to larger attraction and thereby an 
enhanced softening of the EOS. 



We have used the results of Baym, Pethick and Sutherland |32] to describe the crust 
consisting of leptons and nuclei at the low-density (% < 0.001 fm~ 3 ) EOS. For the mid- 



density regime (0.001 < % < 0.08 fm -3 ) the results of Negele and Vautherin []33| are 
employed. Above this density, the EOS for the relativistic models have been adopted. It is 
worth mentioning here that though we have treated K condensation as a second order phase 
transition, for higher U R (n ) values the incompressibility K = 9dP/driB at the threshold is 
found to be negative (see Fig. 5). This represents a first order phase transition where we 
have employed the Maxwell construction to maintain a positive compressibility. 

The static neutron star sequences representing the stellar masses M/M Q and the cor- 
responding central energy density e c are shown in Fig. 6 for nucleons-only matter (solid 
line) and matter with further inclusion of K condensate (dashed line) for several values of 
Ux(no). With the occurrence of K condensation, the maximum masses of neutron stars are 
reduced due to softening of the EOS. The maximum masses M max and the central densities 
Wccnt = n cent /no of the stars are listed in Table III. Interestingly, it is observed that for 
Ur^o) > —140 MeV, the threshold densities for neutral K° condensate occur before the 
central densities for the maximum mass stars. This implies that a significant region of these 
maximum mass stars will contain K° condensate along with K~ condensate. Note that for 
large values of U R (n ), the requirement of Maxwell construction causes a negligible variation 
of the masses (and radii) of the stars with central density. 

Apart from the K optical potential depth, the critical densities for antikaon condensation 
depend sensitively on the nuclear equation of state at high density. The GM1 set is obtained 
from a Walecka-type Lagrangian with a self-interaction term for the scalar meson. Because 
of the linear density-dependence of the vector w-meson field, it was found that the vector 
potential overestimates the relativistic Bruckner-Hartree-Fock (RBHF) results at high den- 
sities. Bodmer [33] first proposed an additional nonlinear tu-meson term in the RMF model 
of the form 

O = lg 4 {u^f . (28) 



Later, it was found by Sugahara and Toki |35| that this modification led to a reasonable 



agreement with the RBHF results. The model parameters obtained [53] by fitting exper- 
imental data for binding energies and charge radii of heavy nuclei is referred to as TM1. 
The parameters of this set are given in Table I. The a — K coupling constants g a x for 
the TM1 set listed in Table II are found to be much smaller than GM1. This stems from 
rather large uj — K potential of Vjf(no) = 91 MeV for the TM1 set compared to 72 MeV 
for GM1. In the TM1 set, for K depths of U R (n ) = -100(-180) MeV, we obtain KN 
average scattering lengths and kaon optical potential depths of o>kn = — 0.829(— 0.373) fm 
and U K (n ) « +82(+3) MeV. In this set even for U R (n Q ) = -180 MeV, the kaon optical 
potential depth is repulsive due to rather large V^f (no) in contrast to GM1 set. 



10 



The equation of state for nucleons-only star matter with and without antikaon condensa- 
tion is shown in Fig. 7 for the TM1 set. In contrast to the GM1 set, the EOS for np-matter 

1 /3 

at high densities is much softer here. This is due to n B variation of the vector potential 
compared to simple n# dependence in the GM1 set. Due to a soft EOS, the meson fields in 
the TM1 set vary rather slowly with density leading to delayed appearance of K condensate 
(see Table III). The stellar sequence for this set is shown in Fig. 8. It is found that the 
masses of the stars with only nucleons and leptons are smaller than those in GM1. On the 
contrary, with K condensate the M max of stars for Ux(n ) > —120 MeV are in fact larger 
in the TM1 model compared to the GM1 set. The delayed formation of K in this model 
suppresses the softening effect of antikaons. In particular, only for K depth as large as —160 
MeV, a K° condensate can be formed in the limiting mass star. 

The above study suggests that a stiffer EOS is more efficient for antikaon condensation 
in neutron star matter. The small values of no = 0.145 fm -3 and m* N /m^ = 0.634 with 
large a asy = 36.9 MeV in the TM1 set should have, in principle, provided a stiffer EOS [|J. 
However, the nonlinear cu-meson term in this model entails a substantial softening of the 
EOS at high densities. To investigate the effects of a rather stiff EOS on K condensation 
with acceptable nuclear matter properties, we have generated another set of parameters 
using the Lagrangian of Eq. (1) (i.e. excluding the nonlinear cu-meson term of Eq. (28)) 
and reproducing the saturation properties of the TM1 set. We call this set as GMT and 
it is presented in Table I. It is to be noted that the coefficients in the sets GM1 and 
GMT are negative. This is associated with the well-known problem that the scalar potential 
is unbounded from below. We have however found that in these effective field theoretical 
models the equations of state for the two sets are continuous (as are exhibited) over a wide 
density range relevant to the neutron star interior. The g a x coupling constants presented 
in Table II for the GMT set are nearly the same as TM1. This arises from similar values of 
uj — K potential at saturation density although they should differ considerably at higher 
densities. For the same reason, the values of K N scattering lengths and K optical potential 
depths at are almost identical in the TM1 and GMT sets. 

In Fig. 9, we display the variation of effective kaon mass m* K (top panel) and K conden- 
sate energies (bottom panel) for various values of U^(n ) for the GMT set. This should be 
contrasted with Figs. 3 and 4 for the GM1 set. It is found that m* K has a smaller decrement 
with density even for Ug(no) = —180 MeV. Therefore, the variation of tu-meson field is 
chiefly responsible for the large drop in uJk-,k° i n the GMT set. The EOS and the stellar 
sequences are shown in Figs. 10 and 11 with and without K condensation. The results for 
the GMT set are summarized in Table IV. The M max for nucleons-only stars is largest in the 
GMT set. Since the EOS is very stiff, the chemical potentials for nucleons increase rapidly 
with density which shift the threshold densities for K condensation to much smaller values. 
We find that K° condensation occurs well inside the maximum mass stars for all values of 
Uk{tiq), except for —100 MeV. To delineate the effects of K° condensation from K~, we also 
give in Table IV the maximum masses and corresponding central densities for stars allowing 
only K~ condensation. It is found that the softening due to K° condensate could reduce the 
maximum mass from that of the K~ condensate further by ~ 17% at U^no) = —180 MeV. 
The present investigation points to the fact that a stiff EOS (as for the GMT set) favors 
strongly the formation of K° condensate in neutron stars. In other words, from G-parity one 
can infer that a soft EOS would enhance kaon production. Indeed, more kaon production 



11 



was found [p6| |38|] in heavy ion collisions for a soft EOS rather than a stiff one. 

We have discussed in Fig. 3 that K~ condensation which occurs at earlier density than 
K° reduces the magnitude of the vector potentials. Consequently, the effective energy u^o 
of K° is enhanced shifting its threshold to a higher density. This effect is illustrated in Table 
V where the stellar properties of stars having only K° condensate is considered. It is evident 
that K° formation takes place much earlier in the absence of K~ condensate (see Tables III 
and IV). Since u cr (K°) < u cent , sizeable amount of K° condensate would be present in stars 
for all K optical potential depths in the GM1 and GMT sets. However, u cr (K°) being higher 
than u CT (K~), and moreover, K° cannot replace leptons in maintaining charge neutrality, 
the EOS is stiffer and the maximum masses of the stars are larger here than those with K~ 
condensate only. 

In Fig. 12, the mass-radius relationship is shown for nucleons-only stars and for stars 
with different values of Uj^(rio) for the GM1, TM1, and GMT sets. It is found that in all 
models the limiting mass stars without K condensate have the largest maximum masses with 
smallest radii. With increasing K optical potential depth although the maximum masses 
decrease, the corresponding radii are however found to be similar for most of the stars in 
the TM1 and GMT sets. In contrast, stars with smaller masses have larger radii in the GM1 
set. 

We now briefly consider the case when hyperons are allowed in addition to nucleons. At 
higher densities when the Fermi energy of nucleons exceeds the effective mass of a hyperon 
minus its associated interaction, the conversion of nucleons to hyperon is energetically fa- 
vorable. The hyperon-meson coupling constants are obtained as discussed before. Unless 
otherwise mentioned we use an attractive E well-depth of —30 MeV. With this choice the A 
and S~ are the first two strange particles that appear at roughly the same density of 2no in 
all the sets. More massive and positively charged particles than these appear at higher den- 
sities. Since the conversion of nucleons to hyperons relieves the Fermi pressure of nucleons, 
the equation of state is softened. The EOSs for star matter with only nucleons and hyperons 
are depicted by dash-dotted lines in Figs. 5 and 7 for the GM1 and TM1 sets. The maximum 
masses and the corresponding central densities are given in Table III. It is evident from the 
figures that the hyperons emerge at densities much earlier than the threshold density for K 
condensate even for U^{no) = —180 MeV. Since we have seen that a soft EOS prevents K 
condensation, the hyperons therefore postpone K formation to higher densities. This effect 
is further accentuated because the negatively charged S~ and S~ can replace electrons in 
maintaining charge neutrality. The resulting drop in /i e would mean a smaller value of ujk- 
is required for K~ condensation which leads to higher threshold density. However, we find 
that in presence of hyperons the rate of decrease of the isovector potential is not so severe 
compared to that of /i e . Therefore, in a hyperonic medium the shift in the threshold density 
for K° condensate is much smaller than that of K~ . For example, in the GM1 set, the 
critical densities are u cr (K~) = 6.00n and u cr (K°) = 6.61n for Ux(n ) = —120 MeV. 
(This is to be compared with the results of Table III for stars with nucleons and antikaons.) 
These values are however beyond the central density of 5.35no for maximum mass star con- 
taining all baryons but no antikaons. On the other hand, for Ux{no) = —140 MeV, the 
critical densities of u cr (K~) = 3.41no and u CI (K°) = 4.42no imply that both K~ and K° 
condensates will be present in hyperonic stars. This means that even though hyperons shift 
the threshold densities of K~ and K° condensates to higher values, the possibility of finding 
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hyperon stars with K° condensate (in association with K ) is enhanced for large K optical 
potential depths. 



The standard mean field model used in the hyperon sector was found [39 to be inad- 
equate to describe the strongly attractive hyperon-hyperon interaction observed in double 
A-hypernuclei. Since the core of a neutron star can be hyperon-rich [Q, the hyperon-hyperon 
interaction is expected to be important. This may be accounted by considering two addi- 
tional mesons, the scalar meson /o(975) (denoted by a* hereafter) and the vector meson 



0(1020) which couple only to the hyperons (Y) [|13| , p9| . The corresponding Lagrangian is 

B 

+\ (d M <7W - m>* 2 ) - \<P^ + \ml<P,r • (29) 
The <p — Y couplings are obtained from the SU(6) relation 

2^a = 2^2 = g^E = -^-g^N • (30) 

The (j* — Y couplings are obtained by fitting them to a well-depth U Y for a Y in a Y -bath 
13| , |3"9"|| . Note that the nucleons do not couple to the strange mesons, i.e. g a *N = g^N = 0. 



The EOS and the star sequence are shown for this situation in Figs. 5 and 6 (dotted lines) 
for the GM1 set without antikaons. The additional attraction imparted by the a* field makes 
the EOS softer at moderately high densities. At very high densities, the EOS merges with 
that without the strange mesons due to dominance of the repulsive <fr field over the a* field. 
The maximum star mass obtained here is 1.721M at a central density of 5n . 

The rather large values of the meson fields (cr*,0) in a hyperonic star should influence 
antikaon condensation densities. From the strangeness content of K (s-quark content) it 
is evident that the a* field is attractive and the field is repulsive for antikaons as also 
for the hyperons. As in Ref. ||13|| , the <p — K coupling is obtained from the SU(3) relation, 
g<t>K = 6.04/v / 2, and the a* — K coupling from the /o(975) decay, g^K = 2.65. With 
these couplings and for U^{no) = —140 MeV, we find that the relatively strong repulsive <fr 
field in a hyperon-rich medium prevents antikaon condensation even at densities as large as 
7.40rio (in the GM1 set) where the effective nucleon mass m* N drops to zero P^|JT3]| . Only 
for large K optical depth of —160 MeV, and thereby large g a K value, critical densities of 
u CT (K~) = 2.76n and u cr (K°) = 4.13n are reached. 

The above discussions of hyperonic star matter are confined to S-depth of —30 MeV. 
Repulsive depth of Ug(n ) ~ +30 MeV has been also estimated [^l]]. On refitting the scalar 
coupling g CT £ to this depth and limiting to the Lagrangian of Eq. (1) only, we find that in 
all the models S~ formation occurs at a much higher density and the resulting EOS is now 
stiffer compared to those with attractive S-depth. For example, in the GM1 set, a star of 
maximum mass 1.798M is formed with a central density of 5.13n . All the S particles do 
not exist in this star while the H~ occurs at a early density of 2.6uq. 

In this paper, we have employed the relativistic mean field model where the nucleon- 
meson coupling constants are fitted to certain bulk properties at normal nuclear matter 
density. The model is then extended to high density regime to investigate the effects of 
antikaon condensation. The underlying assumption in the RMF models is that the meson 
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field operators in the Euler-Lagrange equations are replaced by their mean values because 
the source terms in those equations increase with baryon density. Also, the coupling con- 
stants of the RMF models are density independent. There are other approaches where 
gross properties of dense matter relevant to neutron stars are calculated in the framework of 
nonrelativistic Brueckner-Hartree-Fock and variational chain summation (VCS) calculations 



using modern nucleon-nucleon interaction [ 18, 40, 41 1. Recently, Akmal, Pandharipande and 



Ravenhall (APR) []TH] have studied neutron star properties in the VCS method using a new 
Argonne Vis nucleon-nucleon interaction (which fit all the data of the Nijmegen data base) 
with relativistic boost corrections and a fitted three-nucleon interaction. This calculation 
represents at present perhaps the most sophisticated many-body approach to dense matter. 
The bulk properties obtained in the APR model are binding energy E/B = 16 MeV, nuclear 
symmetry energy a asy = 33.94 MeV, and incompressibility K = 266 MeV at a saturation 
density of n = 0.16 fm -3 . To gauge the uncertainties involved in the extrapolation of the 
RMF model at high density, we compare the nucleons-only matter results with those of the 
APR model. For this purpose we use the model which contains nonlinear interaction for 
the scalar cr-meson only. (Note that the TM1 model which includes also self-interaction of 
tu-meson is fitted to finite nuclear properties.) The five parameters in this RMF model are 
fitted to the above mentioned (four) saturation properties plus the effective nucleon mass 
of m* N /rriN = 0.70 and 0.78 for the two representative cases considered here. In Fig. 13 we 
show the results for pure neutron matter (PNM) and symmetric nuclear matter (SNM) in 
the RMF model (solid lines), and in the APR model with pion condensation (dashed lines) 
and their low density extrapolations without a pion condensed phase (dotted lines) Ul8| . It 
is clear from this figure that the agreement between the two models breaks down at high 
density. The APR model takes into account possibly all leading many-body correlation ef- 



fects. It was also demonstrated [0,0 that the strong tensor correlation gives rise to neutral 
pion condensation in PNM and SNM at densities of 0.20 fm~ 3 and 0.32 fm~ 3 , respectively. 
The RMF model, however, cannot address the question of density dependent correlation 
in dense matter, and moreover the coupling constants are density independent. This may 
be crucial to the deviation of the RMF results from the APR calculations. Consequently, 
the threshold densities of antikaon condensation may be affected. On the other hand, the 
advantage of the RMF models is that relativity is in-built and this causes the EOS to be 
causal. However, in the APR model the EOS becomes superluminal at densities close to the 
maximum masses of the stars. 

We have computed nucleons-only star masses using this parameter set in the RMF model. 
For the case m* N /mN = 0.70, the maximum masses of the neutron star matter (pure neutron 
matter) are M max = 2.29 (2.45)M at central densities of n cent = 0.92 (0.80) fm~ 3 . For the 
case m* N /m^ = 0.78, the respective values are M max = 2.00 (2.2O)M at n cent = 1.09 (0.91) 
fm -3 . Comparing with the respective results of M max = 2.20 (2.21)M in the APR model 
with pion condensation, we find that in the RMF model the maximum mass of a neutron 
star is more sensitive to the amount of proton fraction. It was shown [42| that a neutron star 
may cool rapidly by the so-called direct URCA process in which neutrinos generated in the 
stars interior carry away the energy. The threshold density for this process depends on the 
symmetry energy which in turn is determined by the proton fraction. For the APR model 
the threshold is at a density of % = 0.78 fm~ 3 corresponding to a star of mass ~ 2.OM . In 
the RMF model, for m* N /m N = 0.70 and 0.78, the threshold densities are respectively, 0.59 
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and 0.61 fm -3 , which correspond to neutron stars of masses 2.09 and 1.78M . The relatively 
rapid rise of symmetry energy with density manifests in a smaller threshold density for the 
direct URCA process in the RMF models. 

It may be worth mentioning that the role of nucleon-nucleon and kaon-nucleon corre- 
lations in kaon condensation in neutron star matter was studied by Pandharipande et al. 
||43|| . The kaon energy for low density matter was calculated from Lenz potential while at 
high density the Hartree potential was used. The strong NN and KN correlations in the 
Hartree potential were shown to lead to a dramatic reduction of the antikaon attraction 
at high densities that could raise the threshold density for antikaon condensation to much 
higher values. Recently, kaon energy in neutron matter has been calculated analytically ^3 



by solving the Klein-Gordon equation in the Wigner-Seitz cell approximation. It was found 
that the transition from the low density Lenz potential to a high density Hartree potential 



occurs at ~ 4no- However, in a chiral perturbation theory, Waas et al. [45| found that 



short-range correlations has only a moderate effect on the antikaon condensation densities. 



IV. SUMMARY AND CONCLUSIONS 

We have investigated antikaon, K~- and A^-meson, condensation in neutron stars within 
a relativistic mean field approach where the interaction between the baryons and antikaons 
are generated by the exchange of a, u, and p mesons. Three different parameter sets (GM1, 
TM1, and GMT) have been exploited in this calculation. In the GM1 and GMT sets the 
Lagrangian contain self-interaction for a-meson only, while the set TM1 incorporates also 
nonlinear interaction for a;-meson. The softest nuclear equation of state follows from TM1 
set whereas the stiffest one from GMT set. We find that the critical densities for K~ and in 
particular K° condensations depend sensitively on the choice of antikaon optical potential 
depth and more strongly on the nuclear equation of state. The threshold density of K° 
always lie above that of K~ meson. With the appearance of K~ and K° condensates, 
the overall equation of state becomes softer compared to the situation without antikaon 
condensation. This leads to a reduction in the maximum star masses. With the softest EOS 
(TM1 set), K° condensate may be formed inside maximum mass stars only when the K 
optical depth is quite large. On the other hand for the stiffest nuclear EOS (GMT set), K° 
can occur well inside the maximum mass stars for rather small values of antikaon optical 
potential depth. The appearance of K~ meson leads to a smaller variation of the vector 
fields with density. The decrease of K° energy with density is therefore attenuated resulting 
in its formation at higher density than in the case when K~ meson is prohibited. 

With the onset of K~ condensation (neglecting K° meson), the proton fraction rises dra- 
matically and even crosses the neutron fraction at some critical density. From the threshold 
conditions for K~ and K° formations, i.e. ojk- = p n — P p and uj^o = 0, it is evident that 
the critical point of identical proton and neutron densities (Fermi momenta) is precisely 
the density where K° formation sets in. In earlier model studies of K~ condensation only 
(where K° condensation were ignored), identical values of neutron and proton abundances 
were observed at high densities after K~ condensation. It may therefore be instructive to 
revisit the antikaon condensation scenario including K° meson in these models. With the 
onset of K° condensation, there is a competition in the formation of K~ —p and K° — n pairs 
resulting in a perfectly symmetric matter for nucleons and antikaons inside neutron stars. 
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As the proton fraction is much enhanced in this situation, it may have profound implication 
on the cooling properties of neutron stars. 

With the hyperon-meson couplings constants determined from hypernuclear data and 
SU(6) symmetry relation, the hyperons are formed earlier than antikaons. Since the EOS 
becomes softer, antikaon condensation can only occur at high densities. The negatively 
charged hyperons cause a larger reduction in the electron chemical potential compared to 
that of the isovector potential. This indicates that in contrast to nucleons only stars, K° 
has a higher formation probability in hyperonic stars, but only for relatively larger values 
of K optical potential depth. 

For all the cases studied here, the maximum masses of the stars are found to be larger 
than the precise current observational lower limit of 1.44M imposed by the larger mass of 
the binary pulsar PSR 1913 + 16 | 4"6| . The recent discovery of high-frequency brightness 
oscillations in low-mass X-ray binaries provides a promising new avenue for determining 
masses and radii of neutron stars [|7|. Large neutron star masses with an upper limit of 
M ~ 2.0 — 2.3M have been extracted from kilohertz quasi-periodic-oscillations (QPO). 
The present investigation suggests that indeed both K~ and K° condensates are possible in 
stars of such large masses. On the other hand, hyperon degrees of freedom which lower the 
maximum mass considerably, would be clearly excluded in these massive stars. In contrast, 
also very small masses of M « 1.35 ± O.O4M have been accurately determined in binary 
pulsars |48fl. It is clear that exotic degrees of freedom like antikaons and hyperons must 



occur in these pulsars at few times the normal nuclear matter density. 

It has been inferred that large magnetic field ~ 10 18 G could exist at the core of neutron 
stars [|9|. Though such large interior field is not accessible to direct observation, their 
existence could influence certain gross properties of neutron stars, such as cooling properties 
of the stars [HTJ] . Recently, there has been a considerable debate about whether a charged 



Bose gas in a constant magnetic field could exhibit a Bose-Einstein condensation [51-54 



Some calculations ruled out such a possibility [ ]5"T| , |52"[ whereas others predicted that it could 
occur in a magnetic field []53|,|54]|. The formation of a Bose-Einstein condensate of neutral 



Bose gas such as a system of K in a strong magnetic field may be an interesting possibility. 
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TABLES 



TABLE I. The nucleon-meson coupling constants in relativistic mean field models. The pa- 
rameter set GM1 is taken from Ref. [20] where the baryons interact via nonlinear cr-meson and 
linear w-meson exchanges. The coupling constants are obtained by reproducing the nuclear matter 
properties of binding energy E/B = —16.3 MeV, baryon density no = 0.153 fm~ 3 , asymmetry 
energy coefficient o a sy 

= 32.5 MeV, incompressibility K = 300 MeV, and effective nucleon mass 
m* N /mN = 0.70. The hadronic masses are mjv = 938 MeV, m c = 550 MeV, = 783 MeV, and 
m p = 770 MeV. The parameter set TM1 is obtained from Ref. [35] which has nonlinear exchanges 
in both a and u mesons. The nuclear matter properties in TM1 set are E/B = —16.3 MeV, 
n = 0.145 fur 3 , a asy = 36.9 MeV, K = 281 MeV, and m* N /m N = 0.634. All the hadronic masses 
in this model are same as GM1 except for cr-meson which is m a = 511.198 MeV. The parameter set 
GMT corresponds to nonlinear a and linear w-meson exchanges as in GM1, however, the nuclear 
matter saturation properties and cr-meson mass of the TM1 set are used. All the parameters are 
dimensionless, except gi which is in fm _1 . 

9aN_ 9uiN 9pN_ 02 93 ff4 

GM1 ~ 9.5708 10.5964 8.1957 12.2817 -8.9780 

TM1 10.0289 12.6139 4.6322 -7.2325 0.6183 71.3075 

GMT 9.9400 12.2981 9.2756 10.5745 -24.1907 



TABLE II. The couplings constants for the antikaons K to the cr-meson, g a x, for various values 
of K optical potential depths Ur (in MeV) at the saturation density. The results are for the GM1, 
TM1, and GMT sets. The couplings constants of K to the vector oj and p mesons are fixed by 
SU(6) relations. 



Ur 


-100 


-120 


-140 


-160 


-180 


GM1 


0.9542 


1.6337 


2.3142 


2.9937 


3.6742 


TM1 


0.2537 


0.8384 


1.4241 


2.0098 


2.5955 


GMT 


0.2435 


0.8220 


1.4015 


1.9800 


2.5596 
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TABLE III. The maximum masses M max and the corresponding central densities 
Ucent = n cen t/no for nucleons-only (np) star matter and for stars with further inclusion of hy- 
perons (npH). The results are for the GM1 and TM1 set. For neutron star matter with nucleons 
and antikaons (npK), the critical densities for the K~ and K° condensation, u cr (K~) and u CI (K°), 
and results for M max and ii ccnt at various values of K optical depth Ur (in MeV) at saturation 
density are also given. 



GM1 


TM1 




Ur 


u CT (K ) 


u CT (K°) 


^■ccnt 


M max /M 


u CT (K ) 


u CT (K°) 


^ccnt 


M max /M 


np 








5.63 


2.364 






5.97 


2.179 




-100 


3.45 


5.51 


5.17 


2.211 


4.15 


11.12 


5.67 


2.142 




-120 


3.05 


4.83 


5.19 


2.077 


3.55 


9.13 


5.55 


2.083 


npK 


-140 


2.71 


4.19 


4.75 


1.856 


3.05 


7.43 


5.65 


1.986 




-160 


2.43 


3.59 


3.59 


1.551 


2.67 


5.99 


6.37 


1.857 




-180 


2.19 


3.07 


3.09 


1.217 


2.37 


4.81 


6.01 


1.641 


npH 








5.35 


1.784 






5.87 


1.653 



TABLE IV. Same as in Table III, but for the GMT set. The maximum masses and correspond- 
ing central densities for stars with K~ condensation only without any K° formation are given in 



the parentheses. 

Ur u ct (K-) u cr (K°) 

^■ccnt 

M max /M 

np - - - 4.91 2.667 

-100 3.23 4.95 4.69 2.546 

-120 2.87 4.45 4.77(4.76) 2.444(2.446) 

npK -140 2.57 3.97 4.75(5.17) 2.279(2.309) 

-160 2.31 3.51 4.78(5.52) 2.012(2.162) 

-180 2.09 3.05 4.37(4.67) 1.584(1.929) 
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TABLE V. The critical densities for K° meson condensation, u cr (K°) (where u = Hb/tio), 
for various values of K optical potential depths U k (jiq) (in MeV) in star matter with nucleons 
and K° mesons only, neglecting K ~ formation. The maximum masses of the stars M max and the 
corresponding central densities u cen t are given. The results are for the GM1, TM1, and GMT 
parameter sets. 



GM1 


TM1 


GMT 


Uk 


u CI (K°) 


^ccnt 


M max /M 


u cr (K°) 


^ccnt 


M max /M Q 


u CT (K°) 


^cent 


M max /M 


-100 


4.80 


5.57 


2.349 


9.35 


5.97 


2.179 


4.32 


4.94 


2.655 


-120 


4.27 


5.35 


2.303 


7.70 


5.97 


2.179 


3.91 


4.80 


2.623 


-140 


3.76 


5.02 


2.205 


6.32 


5.97 


2.179 


3.48 


4.64 


2.553 


-160 


3.32 


4.85 


2.031 


5.22 


6.04 


2.170 


3.13 


4.41 


2.423 


-180 


2.95 


5.35 


1.700 


4.33 


5.73 


2.105 


2.79 


4.08 


2.194 
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FIG. 1. The mean meson potentials for nucleons and the electrochemical potential versus the 
baryon density n^/no in the GM1 set for nucleons-only star matter (solid lines) and for matter with 
further inclusion of antikaon, K~ and K°, condensation (dashed lines). The K optical potential 
depth is Uk = —120 MeV at the normal nuclear matter density of no = 0.153 fm~ 3 . Above the 
critical density for K° condensation, n cr (K°) > 4.83no, the electrochemical potential fi e coincides 
with -V*? = -g p NP03- 
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n B /n 

FIG. 2. The proper number densities rii of various compositions in neutron star matter in the 
GM1 model. The results are for nucleons-only matter (top panel), matter with further inclusion of 
K~ condensation (central panel), and matter where also K° condensation is considered (bottom 
panel). The K optical potential at normal nuclear matter density is = —120 MeV. 
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FIG. 3. The variation of effective antikaon mass m* K /rriK as a function of baryon density n^/no 
for star matter with nucleons and antikaons in the GM1 set. The different curves from the top to 
bottom on the left side of the graph correspond to K optical potential depths at normal nuclear 
matter density of Ug = -100, -120, -140, -160, -180 MeV. 
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n B /n 



FIG. 4. The effective energy of K (solid lines) and K° (dashed lines) versus baryon density in 
the GM1 set. The different curves for each antikaon have the same meaning as in Fig. 3. 
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FIG. 5. The equation of state, pressure P vs. energy density e in the GM1 set. The results 
are for nucleons-only (np) star matter (solid line) and for matter with further inclusion of K~ 
and K° condensation (dashed lines) for antikaon optical potential depths at normal density of 
lift = —100, —120, —140, —160 MeV. The equations of state for stars with nucleons and hyperons 
without any antikaons (npH) are shown for Lagrangian of Eq. (1) (dash-dotted line) and with 
further inclusion of Eq. (29) (dotted line). 
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FIG. 6. The neutron star sequences near the limiting mass in the GM1 set. The different curves 
have the same meaning as in Fig 5. The filled circles correspond to the maximum masses, and the 
arrows indicate the minimum mass stars that possess a K° condensate at their centers. 
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FIG. 8. Same as Fig. 6 but for the TM1 set. The antikaon optical potential depths are same as 
in Fig. 7. 
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FIG. 9. In the top panel the effective antikaon mass m* K /niK, and in the bottom panel the 
effective energy of K~ (solid lines) and K° (dashed lines) versus baryon density in the GMT set. 
The different curves from the top to bottom on the left side of the graph in each panel correspond 
to U R (n ) = -100, -120, -140, -160, -180 MeV. 
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FIG. 10. Same as Fig. 5 but for the GMT set. The values of U^(no) are same as in Fig. 9. 
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FIG. 11. Same as Fig. 6 but for the GMT set. The values of U^(no) are same as in Fig. 9. 
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FIG. 12. The mass-radius relation for neutron star sequences for nucleons-only stars and for stars 
with further inclusion of K~ and K° condensation for the GM1 (top panel), TM1 (central panel), 
and GMT (bottom panel) parameter sets. The filled circles correspond to the maximum masses of 
the stars without and with K condensate. For the latter case, the circles from right to left in each 
panel are obtained with U R (n ) = -100, -120, -140, -160, -180 MeV. 
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FIG. 13. The pure neutron matter (PNM) and symmetric nuclear matter (SNM) energies as a 
function of baryon density. The calculations are in the RMF model (solid lines) with effective 
nucleon mass of m* N /mN = 0.70 and 0.78, and in the APR model [18] with pion condensation 
(dashed lines) and their low density extrapolations without a pion condensed phase (dotted lines); 
for details see text. 
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